%------------------Klauder子波的时域表达式---------------%
function w = Klauder(dt,N,T,f1,f2,fa)
    t=-N/2*dt:dt:N/2*dt;
    k=(f2-f1)/T;
    p=0.5*(f2+f1);
    s1=sinc(k*(T-t).*t).*(T-t).*cos(2*pi*p*t);
    s2=hilbert(s1);
    w=real(s2)*cos(fa)+imag(s2)*sin(fa);
end